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Abstract 



< 

\ Bosonization of the two-dimensional QCD in the large Nc limit is performed in the 

framework of Hamiltonian approach in the Coulomb gauge. The generalized Bogoli- 
ubov transformation is applied to diagonalize the Hamiltonian in the bosonic sector 
\ of the theory, and the composite operators creating/annihilating bosons are obtained 

■ in terms of dressed quark operators. The bound state equation is reconstructed as 

a result of the generalized Bogoliubov transformation, and the form of its massless 
Q \ solution, chiral pion, is found explicitly. Chiral properties of the theory are discussed. 

o\ 

Two-dimensional quantum chromodynamics {QCD2) in the large Nc limit was first con- 
sidered many years ago but still remains popular in studies of various aspects of strong 
interactions. The reason for this is two-fold. There are no transverse gluons in two di- 
mensions, so the theory is relatively simple; moreover, in the large Nq limit only planar 
^ • graphs are to be summed, and a simple diagrammatic approach can be developed. On the 
other hand, this theory does have nontrivial content exhibiting both confinement and chiral 
^ ! symmetry breaking. 

^ I Most of the studies in QCD2 were performed in the light-cone gauge, which consider- 

ably simplifies the spectrum calculations but yields the perturbative vacuum. Alternative 
approach based on the Coulomb gauge = is more technically involved 0, and it is not 
a straightforward exercise to demonstrate the equivalence of the two formulations. In par- 
ticular, it appears that the vacuum is nontrivial in the Coulomb gauge version, and nonzero 
quark condensate exists for massless quarks The latter feature is confirmed by the sum 
rules approach in the light-cone gauge. The confused situation was resolved to large 
extent in the formulation on finite intervals |^, and the equivalence was demonstrated. 

In the present paper we study the Hamiltonian formulation of QCD2 in the Coulomb 
gauge. In contrast to |Q the theory is quantized at the ordinary time hypersurface. Our 
main purpose is to investigate some peculiar properties of mesonic wave functions. It will 
be shown that the Bogoliubov- Valatin approach offers the most natural setting for such 
studies. Finally we perform the generalized Bogoliubov transformation and reformulate the 
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theory in terms of the effective mesonic degrees of freedom. The presented approach allows 
straightforward calculations of any matrix elements of quark operators, and we demonstrate 
this property at the example of evaluation of the pionic decay constant. 
The Lagrangian of QCD2 has the form 

L{x) = -\f;Ax)F;,{x) + q{x){ib - m)q{x) (1) 

where D = {d^—igAp;°')'y^, and our convention for 7 matrices is 70 = 0-3, 71 = i(T2, 75 = 7o7i- 
The large Nc limit implies that g'^Nc remains finite. 

Gluonic propagator in the Coulomb gauge = takes the form Do(){ko,k) = —1/k'^, 
and the infrared singularity is regularized by the principle value prescription yielding the 
linear confinement 

D^'oi^o -yo,x-y) = -'-6'''\x - y\6ixo - y^). (2) 
Hamiltonian can be obtained by standard methods and reads 

H = j dxq+ix) (^-^7'^ + ^^7°) q{x) - y / dxdy q^{x)ty{x)q'^{y)t''q{y)^-^^, (3) 
where the quark and antiquark fields are defined by 

qi{xo, = J ^ [u{k)bi{xo, k) + v{-k)d^{xo, -k)] e'^'' (4) 



6,(A;)|0) = rf,(-A;)|0) =0 (5) 

with 

u{k)=T{k)Q vik)=Tik)Q (6) 

Parameter 6{k) has the meaning of Bogoliubov-Valatin angle describing rotation from 
bare to dressed quarks. Following [0] Hamiltonian (§) can be normally ordered in the basis 
of fermionic operators (^: 

H = LNc£,+ : H2: + : H^:, (7) 
where is the vacuum energy density {L being the one-dimensional volume of the space) 

S. = J ^Tr{{^'p + m7°) A-ip)} + ^^^r {A4k)A4p)} (8) 



9' f.r M 



4tt \ NcJ Nc-^oo 4n 



c 



Part : H2 : is bilinear in quark fields: 

d 



: H2:= j dx: qfi^x) i "^7^^ + "^7° j <li{x) ■ 



1 j dxdy^-^^ J dk : qtix) [A^{k) - A^k)] q,{y) : e^'^^-"^) (9) 



and the projection operators A-t are introduced as 

A± = T{k) 



(10) 



The angle 9{k) is defined from the condition of diagonahzing the bihnear in quark fields 
part of the Hamiltonian : H2 :. This condition gives a system of integral equation for 9{p) 
and quark dispersion E{j)) 



E(p)cos9(p) = m + -f - — cosOik) 
^ J (p — k)^ 



E{p)sin9{p) = P + ^ j~ 



dk 



{p - k) 



;sin9{k), 



which leads to a gap equation for the angle 9{p) 



p cos 9{p) — m sin 9{j)) 
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dk 



{p — k)' 



sm[9{p) - 9ik)]. 



(12) 



Note that once the "dressed" quarks @, define excitations over the true quark vacuum 
|0), then gap equation ([T^) could be reconstructed as the extremum condition 



69{p) 



0, 



(13) 



which ensures that the vacuum energy is minimal. 

Thus the diagonalized Hamiltonian : H2 : takes the form 



: Ho 



dk 
2^ 



E{k) {bt{k)bi{k) + dt{-k)di{-k)} 



(14) 

-k), E{k) = 



The obvious properties of all solutions of the system (|TT|) are 9{k) = — 
E{—k), and 9{k) 7i/2 when k 00. 

Alternatively equation (|T^ could be found as the solution of the Schwinger-Dyson equa- 
tions for the quark propagator: 



S{po,p) 



(15) 



dkodk 



p — m — Ti{p) 

( M2^°'^(^0' ^)^° = [^(P) cos^(P) - + 7'[^(P) sm^(p) - p]. (16) 

[P — K) 

Gap equation (^) contains all information about one-fermion sector of the theory, but 
not about the interaction between fermions, because the : H^^ : part of the Hamiltonian was 
not involved so far. To proceed further we introduce, as it was done in ||^, the colour singlet 
bilinear operators 



B{p.p') = -^btipMp') D{p,p') = -^dt{-p)d,{-p') 
JNc JNc 



M{p,p') 



^d,{-pMp') M+{p,p') = -^bt{p')dti-p) 
INc JNc 



(17) 
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with commutation relations 

[M(p, p')M-^{q, q')] = {D{q, p)5{p' - q') + B{q\ p')5{p - q)} + {2nyS{p' - q')5{p - q) 



[B{p,p')B{q,q')] 



[D{p,p')D{q,q')] 



2n 



2n 



i2ny6{p'-q')6{p~q), 

iB{p, q')5{p' - g) - B{q,p')5{p - g')) 



(Dip, q')5{p' -q)- Diq,p')6ip - q')) 

Nc Nc^oo 



[B{p,p')M{q,q')] = -^M{q,p')5{p - q') ^ 



[B{p,p')M+{q,q')] 
[D{p,p')M{q,q')] 
[Dip,p')M+{q,q')] 



2tx 



M+{p,q)6ip'-q') ^ 
Nr Nc ^oc 



2tx 



Mip',q')6ip-q) 



27T 



M+{p,q')8{p'-q) ^ 0. 



In terms of these operators Hamiltonian @ can be presented as 

dk 



H = LNcS, + ^Nc j —E{k){B{k, k) + D{k, k)} 

7 [ dpdkdQ r e(p)~e(k) . e(Q - p) - e(Q - k) 

/ 7 TTTT TTTT 2 COS Sm 

2 J {2Ti f{p-kY [ 2 2 

X {m+(p,p - Q)D{k -Q,k) + M+{p - Q,p)B{k - Q, k) 

-Bip,p- Q)M{k -Q,k)- D{p,p- Q)M{k, k - Q)} 

e{p) - e{k) e{Q~p)-9{Q-k) 

+ cos COS 

2 2 
X {b{p - Q,p)B{k, k~Q) + Dip,p- Q)D{k - Q, k) 

+ M+(p - Q,p)M{k -Q,k) + M+{p,p- Q)M{k, k - Q)} 

. e{p) - e{k) . e{Q-p)-e{Q~k) 

+ sm sm 

2 2 

X {b{p,p- Q)D{k, k-Q)+B{p- Q,p)D{k - Q, k) 

+M{p,p- Q)M{k -Q,k) + M+{p - Q,p)M+{k, k - Q)} 



(19) 



(20) 



It can be easily checked that the anzatz ^ 

B{p,p') = ^M^q",p)Miq",p' 



D{p,p') 



iNc 
1 



^ J 27r 



(21) 



M+(p, q")M{p',q" 



satisfies the commutation relations ([T9|), so that in the leading in Nc order Hamiltonian ( pOD 
can be rewritten as 



H = LNcEv + 
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dQdp \ 
(27r)2 

dk 



{E{p) + E{Q - p))M^{p - Q,p)M{p - Q,p) 



{2C{p, k, Q)M+{p - Q, p)M{k - Q, k) 



where 



2 J {p-kf 

+S{p, k, Q) (m{p,p- Q)M{k -Q,k) + M+{p,p - Q)M+{k - Q, k))} 

r( ^(P) - ^(fe) e{Q-p)-9{Q-k) 
6 (p, /c, Q) = cos cos 



Q( i-m ■ W-Ojk) . 9{Q-p)-9{Q-k) 
S [p, k, Q) = sm z sm 



(22) 



(23) 



2 2 

We are now to perform the second Bogoliubov transformation which should diagonalize 
this Hamiltonian. Let us define the new operators as 



(Q) = J^ {M(g - g, q)^l{q, Q) + M+(g, q - Q)^l{q, Q)} 



(24) 



where functions ip^ and parametrize the given transformation and obey Bogoliubov-type 
normalization and completeness conditions (m, n > 0): 



{vl{p,Q)v>7iP^Q)-v>-iP,Q)v-iP,Q)) = 5n 



dp 
2^ 

^ U{p, Q)v^{p, Q) - ^^(j9, g)/;(p, Q) 



5:(y.!;(p,Q)^!^(A;,Q)-y,!!(p,Q)¥,!!(fc,Q)) = 2T,5{p-k) 

n=0 
oo 

$:(y.!^(p,g)^!!(fc,g)-y.'^(p,gX(A;,g)) = o. 

n=0 



(25) 



(26) 



It is easy to see using these conditions that in the Nq oo limit the new operators 
m{Q) and m'^{Q) obey standard bosonic commutation relations 



^n(g)m+ (g')l = 27r5(g - g') s^m, K(g)m^(g')] = U+(g)m+ (g')l = o (27) 
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Straightforward but tedious calculation shows that the transformation (^) diagonalizes 
Hamiltonian (|22|) if functions and y?" are the solutions to the system of equations 

[E{p)+E{Q-p)-Q,]^^{p,Q) 

/d h" 
j-—^ [C{p, k, Q)^+{k, Q) - k, Q)^.{k, Q)] 

[E{p) + E{Q-p) + Qo]^.{P,Q) 

/d h" 
j^-^ [C{p, k, Q)^.{k, Q) - S{p, k, Q)Mk, Q)] , 

^ (28) 
and the resulting Hamiltonian takes the form 

1 +00 „ 

H = LNcS'^ + 7^T. irQniQ) {m+(Q)m„(Q) + m^{Q)m+{Q)} , (29) 

where Q°(Q) is the ra-th positive eigenvalue of the system (pS]). Hamiltonian (|29|) together 
with boson operators (20) comprises our procedure of bosonization. Note that here the vac- 
uum energy density S!^ contains extra contribution which comes from the mesonic operators 
ordering which is of order 0(-^). 

Equations (|28| ) are nothing but those obtained by I.Bars and M.B.Green in for the 
Bethe-Salpeter wave function 

<l>{p, Q) = T{p) f l±^75y,^.(p, Q) + 1^^5y,_(p, Q)) T+{Q - p). (30) 



There is a very important point concerning Bars-Green equations (pSl) . While equations 



(p5D, ( [26| ) are quite natural conditions imposed at the parameters of Bogoliubov-Valatin 
transformations, it is apparently unacceptable to assume such normalization and complete- 
ness conditions for the solutions of bound state equations. In fact the problem is rooted in 
the properties of Bars-Green equations ( pS]) . It can be easily checked that if system ( pSj) is 
rewritten in the form of matrix integral Shrodinger-like equation 

«(:':) ^^-g) (31) 

then the integral operator K is not Hermitian. Fortunately this does not cause a disaster as 
the eigenvalues are real: integrating both sides of ( P^ over p, doing the same for complex 
conjugated equations and taking the appropriate linear combination one arrives at 

E {Ql-Qt) / /(<(p,QK*(p,Q)-^-(p,QVr(p,Q)) =0, (32) 

that gives 

Ql = QT (33) 

together with the orthonormality condition 



^(^!^(p,Q)^!;^(p,Q)-^!!(p,Q)^!^(p,Q)) = 5r, 



dp 
2^ 



(34) 



[vl{p,QW:i{p,Q)-^'l{p,Q)v"l{p,Q)) = 0. 
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Solutions of the system ( pq ) come in pairs: for each eigenvalue Qq with eigenfunction 
(V5+,V^") there exists another eigenvalue —Qq with eigenfunction (v?",'/'"). With this sym- 
metry equation (|3^ can be rewritten in the form (^) where only positive eigenvalues enter. 
Similarly in attempts to construct the Green function for the system (pSD the completeness 
( p6D can be derived. 

From the point of view of Bethe-Salpeter equation the reason for such an unusual norm 
is the following. The matrix equation for the function $ contains projectors (p!OD, so that $ 
is subject to the constraint 

A+ipMp, g)A+(Q -p) = A-ipMp, Q)A-iQ - p) = 0, (35) 



as it is clear from (|30|). Thus both the norm and the completeness are defined in the truncated 



space (|35|) 



It is well known that Bogoliubov-Valatin transformation not only changes the operators, 
the ground state is transformed too. Indeed the boson annihilation operator m„(Q) does not 
annihilate the vacuum |0) defined by equation (^. The explicit expression relating bosonic 
(|f2)) and fermionic (|0)) vacua is rather complicated, but luckily in the large Nq limit the 
matrix elements of quark bilinears calculated in old and new vacua coincide. For example, 
the chiral condensate 

(qq) = {Q\q,{x)q\x)\Q) ^-^^{0\q,ix)q\x)\0) 

= Nc I — rr|7°A_(A;)| = -— / dkcose{k) (36) 

— CO 

is the one found in 

Form ( plf ) for the operator m^{Q) suggests the obvious particle-hole interpretation: the 
wave function of a meson moving forward in time contains two pieces, of a quark-antiquark 
pair moving forward in time with the amplitude ip^ and a pair moving backward with the 
amplitude ip-. No effect of such a kind could ever happen in quark potential models. 

What can be said about the relative size of the two amplitudes ip+ and System 
( pSD was solved numerically in 0, and it was shown that the ip- component is small i) if 
the quark mass is large and ii) for higher excited states. In these two cases quark potential 
model with local linear confinement serves as a good approximation. 

Besides component dies out with the increase of total mesonic momentum Q: in 
the infinite momentum frame {Q oo) it is zero whereas the equation for ip>^ goes to the 
't Hooft equation [jl| after appropriate rescaling, as it was shown in 0. 

In conclusion, let us consider briefly the state which suffers from the effect of backward 
motion in the most dramatic way. It is two-dimensional pion. It was shown in that the 
gap equation had a nontrivial solution in the chiral limit m = 0, and the chiral condensate 
( p6D did not vanish with this solution. So the Goldstone mode should exist in the spectrum. 
Indeed, the set of functions 

</?± {p,Q)=N^ cos ± sm (37) 
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is a solution of the system ( |28D for Qq = y/Q^ and Q > (for Q < ip± <^=^ Vt^-)- Here A^^ 
is the pion norm defined according to the general form (^4|) 



'-00 27T 

or with the solution ( P7| ) substituted one arrives at 

N^iQ) = r — [sin^(p) +sin^(g -p)] = -Q. (39) 
j-oo vr TT 

For Q = ip^{p,0) = (y9_(p, 0) ~ cos^(p), pion spends half of time in the backward 
motion of the pair, and as follows from equation ( |5P| ) such a function has zero norm, as 
it should be for the massless particle at rest. In the opposite limiting case Q ^ oo the 
backward motion part dies out so that 




^lip,Q)-^ J^, 0<P<Q (40) 

coinciding with the Goldstone mode of the 't Hooft equation. This does not mean, neverthe- 
less, that the pionic physics is exhausted by the simple picture of linear confinement in the 
infinite momentum frame. All nontrivial content of the wave function (|37D is concentrated 
in the boundary regions x ^ and x ^ 1, x = p/Q. The same is true of course for the 
QCD2 quantized at the light-cone 0, where quantities like chiral condensate do not come 
out trivially, and one is forced to employ the sum rules approach to arrive at a reliable 
answer. 

With the given Hamiltonian approach it appears straightforward to calculate any ma- 
trix elements of currents between mesonic states. For example, to evaluate the pion decay 
constant defined in the standard way 

{^\4{x)\n{Q)) = UQ,^^, (41) 

we express the axial- vector current J^(a;) = il){x)'^^'^^il){x) in terms of mesonic creation- 
annihilation operators m+, m„ that allows to calculate the matrix element on the l.h.s. 
explicitly and gives for f^, 




U = \—- (42) 



It is instructive to note that for any mesonic state M.n{Q) the analogous matrix element 
{yL \ J^{x) \ Ain{Q)) contains this meson wave function integrated with the pionic one, that 
obviously vanishes for any mesonic state but the pion when the pion norm (^Uj) appears. Thus 
the matrix element (^Tj) is the only nontrivial one and the corresponding decay constants for 
higher excited mesonic states vanish. 

Now we can slightly relax the chiral limit and find the pion mass in the limit of small 
but not exactly vanishing quark mass m. To this end we rewrite the system (^) as a single 
equation for the matrix function $(p, Q) in the form 

Qo*(p, Q) = {l5P + lom)<l>{p, Q) - $(p, Q)(75(<5 - P) + lorn) 
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/fi h 
j——y^ {A+{kMp, Q)A„(g - A:) - A+{pMk, Q)A_(Q - p) (43) 

-A_(A;)$(p, Q)A+(Q - fc) + A_(p)$(fc, Q)A+(Q - p)} . 

On multiplying equation (|i3|) by 7075, taking trace over spinor indices and integrating 
over momentum p one can arrive at 

Qo J ^5j9[7o75$b, Q)]-QJ ^SpI^oHp, Q)] = -2m J ^Sp[^,^p, Q)]. (44) 

With the pion wave function (|37| ) substituted and using definition (pOD the latter equation 
gives a well-known relation 

/^M^ = -2m(gg), (45) 

so that one has for the pion mass 

Ml = 2m / dpcos6{p) ~ rriyfy, (46) 

JO 

which vanishes in the exact chiral limit. 

The last concluding remark concerning Hamiltonian ( p9D is in order. As it could be 
anticipated from the very beginning this Hamiltonian describes free non-interacting mesons, 
whereas the interaction suppressed by powers of A^^ is hidden in the terms present in ( ppj ) but 
omitted in (|29|) . These terms define three and four meson vertices so that their restoration 
gives quite a natural way of considering strong mesons decays and scattering amplitudes. 
This investigation is of a special interest in view of the fact that the pionic wave function 
is found explicitly and hence the "mysterious" infiuence of the qq pair backward motion in 
time upon an exited meson decay into pions can be easily resolved. This work is in progress 
and will be reported elsewhere in the future. 
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